Comparison and optimization of protocols and whole-genome capture conditions for ancient DNA samples.

Ancient DNA (aDNA) obtained from human remains is typically fragmented and present in relatively low amounts. Here we investigate a set of optimal methods for producing aDNA data by comparing silica-based DNA extraction and aDNA library preparation protocols. We also test the efficiency of whole-genome enrichment (WGC) on ancient human samples by modifying a number of parameter combinations. We find that the Dabney extraction protocol performs significantly better than alternatives. We further observed a positive trend with the BEST library protocol indicating lower clonality. Notably, our results suggest that WGC is effective at retrieving endogenous DNA, particularly from poorly-preserved human samples, by increasing human endogenous proportions by 5x. Thus, aDNA studies will be most likely to benefit from our results.

The time and place of origin of South Caucasian languages: insights into past human societies, ecosystems and human population genetics.

This study re-examines the linguistic phylogeny of the South Caucasian linguistic family (aka the Kartvelian linguistic family) and attempts to identify its Urheimat. We apply Bayesian phylogenetics to infer a dated phylogeny of the South Caucasian languages. We infer the Urheimat and the reasons for the split of the Kartvelian languages by taking into consideration (1) the past distribution ranges of wildlife elements whose names can be traced back to proto-Kartvelian roots, (2) the distribution ranges of past cultures and (3) the genetic variations of past and extant human populations. Our best-fit Bayesian phylogenetic model is in agreement with the widely accepted topology suggested by previous studies. However, in contrast to these studies, our model suggests earlier mean split dates, according to which the divergence between Svan and Karto-Zan occurred in the early Copper Age, while Georgian and Zan diverged in the early Iron Age. The split of Zan into Megrelian and Laz is widely attributed to the spread of Georgian and/or Georgian speakers in the seventh-eighth centuries CE. Our analyses place the Kartvelian Urheimat in an area that largely intersects the Colchis glacial refugium in the South Caucasus. The divergence of Kartvelian languages is strongly associated with differences in the rate of technological expansions in relation to landscape heterogeneity, as well as the emergence of state-run communities. Neolithic societies could not colonize dense forests, whereas Copper Age societies made limited progress in this regard, but not to the same degree of success achieved by Bronze and Iron Age societies. The paper also discusses the importance of glacial refugia in laying the foundation for linguistic families and where Indo-European languages might have originated.

Single hybrid population but multiple parental individuals at the origin of parthenogenetic rock lizards Darevskia sapphirina and D. bendimahiensis Schmidtler, & Eiselt Darevsky (1994) endemic to the area of Lake Van in East Turkey.

Among vertebrates, obligate parthenogenesis is only found in Squamata, where it always has a hybrid origin and a few lizard genera contain most of the known hybridogenous parthenogenetic taxa. Parthenogenesis thus seems to be pre-conditioned at the genus level, but it is not clear how often the encounter between two parental sexually reproducing species can result in the parthenogenetic offspring, nor whether the success of such hybridization event requires certain conditions or the specific time frame. To address this question, we studied the rock lizards of genus Darevskia, where a pair of parental species, D. valentini and D. raddei, as well as their parthenogenetic daughter species D. bendimahiensis and D. sapphirina, are found in close proximity NE of the Lake Van in East Anatolia. Using ddRAD-seq genotyping on 19 parental and 18 hybrid individuals, we found that (i) all parthenogenetic individuals from both D. bendimahiensis and D. sapphirina have a monophyletic origin tracing back to a single initial hybrid population, but their current genetic variation is geographically structured; (ii) unlike the most probable paternal ancestor, the genetically closest extant population of the maternal ancestor is not the geographically nearest one; and (iii) in the parthenogens, about 1% of loci carry multiple haplotypes, frequently differentiated by multiple substitutions. This pattern, in addition to biases in the relative frequency of haplotypes of maternal and paternal origin, does not appear compatible with a scenario of the entire parthenogenic clonal population having descended from a single pair of parental individuals. Instead, the data suggest that multiple parental individual ancestries still persist in the parthenogenetic gene pool. This supports the notion that although hybridization leading to parthenogenesis is generally rare at the level of species, it may be more common at the individual/population level once the right conditions are met.

Mobility and kinship in the world's first village societies.

Around 10,000 y ago in southwest Asia, the cessation of a mobile lifestyle and the emergence of the first village communities during the Neolithic marked a fundamental change in human history. The first communities were small (tens to hundreds of individuals) but remained semisedentary. So-called megasites appeared soon after, occupied by thousands of more sedentary inhabitants. Accompanying this shift, the material culture and ancient ecological data indicate profound changes in economic and social behavior. A shift from residential to logistical mobility and increasing population size are clear and can be explained by either changes in fertility and/or aggregation of local groups. However, as sedentism increased, small early communities likely risked inbreeding without maintaining or establishing exogamous relationships typical of hunter-gatherers. Megasites, where large populations would have made endogamy sustainable, could have avoided this risk. To examine the role of kinship practices in the rise of megasites, we measured strontium and oxygen isotopes in tooth enamel from 99 individuals buried at Pinarbasi, Boncuklu, and Çatalhöyük (Turkey) over 7,000 y. These sites are geographically proximate and, critically, span both early sedentary behaviors (Pinarbasi and Boncuklu) and the rise of a local megasite (Çatalhöyük). Our data are consistent with the presence of only local individuals at Pinarbasi and Boncuklu, whereas at Çatalhöyük, several nonlocals are present. The Çatalhöyük data stand in contrast to other megasites where bioarchaeological evidence has pointed to strict endogamy. These different kinship behaviors suggest that megasites may have arisen by employing unique, community-specific kinship practices.

Somatic copy number variant load in neurons of healthy controls and Alzheimer's disease patients.

The possible role of somatic copy number variations (CNVs) in Alzheimer's disease (AD) aetiology has been controversial. Although cytogenetic studies suggested increased CNV loads in AD brains, a recent single-cell whole-genome sequencing (scWGS) experiment, studying frontal cortex brain samples, found no such evidence. Here we readdressed this issue using low-coverage scWGS on pyramidal neurons dissected via both laser capture microdissection (LCM) and fluorescence activated cell sorting (FACS) across five brain regions: entorhinal cortex, temporal cortex, hippocampal CA1, hippocampal CA3, and the cerebellum. Among reliably detected somatic CNVs identified in 1301 cells obtained from the brains of 13 AD patients and 7 healthy controls, deletions were more frequent compared to duplications. Interestingly, we observed slightly higher frequencies of CNV events in cells from AD compared to similar numbers of cells from controls (4.1% vs. 1.4%, or 0.9% vs. 0.7%, using different filtering approaches), although the differences were not statistically significant. On the technical aspects, we observed that LCM-isolated cells show higher within-cell read depth variation compared to cells isolated with FACS. To reduce within-cell read depth variation, we proposed a principal component analysis-based denoising approach that significantly improves signal-to-noise ratios. Lastly, we showed that LCM-isolated neurons in AD harbour slightly more read depth variability than neurons of controls, which might be related to the reported hyperploid profiles of some AD-affected neurons.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data.

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

A genomic snapshot of demographic and cultural dynamism in Upper Mesopotamia during the Neolithic Transition.

Upper Mesopotamia played a key role in the Neolithic Transition in Southwest Asia through marked innovations in symbolism, technology, and diet. We present 13 ancient genomes (c. 8500 to 7500 cal BCE) from Pre-Pottery Neolithic Çayönü in the Tigris basin together with bioarchaeological and material culture data. Our findings reveal that Çayönü was a genetically diverse population, carrying mixed ancestry from western and eastern Fertile Crescent, and that the community received immigrants. Our results further suggest that the community was organized along biological family lines. We document bodily interventions such as head shaping and cauterization among the individuals examined, reflecting Çayönü's cultural ingenuity. Last, we identify Upper Mesopotamia as the likely source of eastern gene flow into Neolithic Anatolia, in line with material culture evidence. We hypothesize that Upper Mesopotamia's cultural dynamism during the Neolithic Transition was the product not only of its fertile lands but also of its interregional demographic connections.

Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity.

Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver, and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence (DiCo). We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell-type composition and also by cell-autonomous expression changes within particular cell types.

MTaxi: A comparative tool for taxon identification of ultra low coverage ancient genomes.

A major challenge in zooarchaeology is to morphologically distinguish closely related species' remains, especially using small bone fragments. Shotgun sequencing aDNA from archeological remains and comparative alignment to the candidate species' reference genomes will only apply when reference nuclear genomes of comparable quality are available, and may still fail when coverages are low. Here, we propose an alternative method, MTaxi, that uses highly accessible mitochondrial DNA (mtDNA) to distinguish between pairs of closely related species from ancient DNA sequences. MTaxi utilises mtDNA transversion-type substitutions between pairs of candidate species, assigns reads to either species, and performs a binomial test to determine the sample taxon. We tested MTaxi on sheep/goat and horse/donkey data, between which zooarchaeological classification can be challenging in ways that epitomise our case. The method performed efficiently on simulated ancient genomes down to 0.3x mitochondrial coverage for both sheep/goat and horse/donkey, with no false positives. Trials on n=18 ancient sheep/goat samples and n=10 horse/donkey samples of known species identity also yielded 100% accuracy. Overall, MTaxi provides a straightforward approach to classify closely related species that are difficult to distinguish through zooarchaeological methods using low coverage aDNA data, especially when similar quality reference genomes are unavailable. MTaxi is freely available at https://github.com/goztag/MTaxi.

Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication.

Sheep were among the first domesticated animals, but their demographic history is little understood. Here we analyzed nuclear polymorphism and mitochondrial data (mtDNA) from ancient central and west Anatolian sheep dating from Epipaleolithic to late Neolithic, comparatively with modern-day breeds and central Asian Neolithic/Bronze Age sheep (OBI). Analyzing ancient nuclear data, we found that Anatolian Neolithic sheep (ANS) are genetically closest to present-day European breeds relative to Asian breeds, a conclusion supported by mtDNA haplogroup frequencies. In contrast, OBI showed higher genetic affinity to present-day Asian breeds. These results suggest that the east-west genetic structure observed in present-day breeds had already emerged by 6000 BCE, hinting at multiple sheep domestication episodes or early wild introgression in southwest Asia. Furthermore, we found that ANS are genetically distinct from all modern breeds. Our results suggest that European and Anatolian domestic sheep gene pools have been strongly remolded since the Neolithic.

Human inbreeding has decreased in time through the Holocene.

The history of human inbreeding is controversial.1 In particular, how the development of sedentary and/or agricultural societies may have influenced overall inbreeding levels, relative to those of hunter-gatherer communities, is unclear.2-5 Here, we present an approach for reliable estimation of runs of homozygosity (ROHs) in genomes with ≥3× mean sequence coverage across >1 million SNPs and apply this to 411 ancient Eurasian genomes from the last 15,000 years.5-34 We show that the frequency of inbreeding, as measured by ROHs, has decreased over time. The strongest effect is associated with the Neolithic transition, but the trend has since continued, indicating a population size effect on inbreeding prevalence. We further show that most inbreeding in our historical sample can be attributed to small population size instead of consanguinity. Cases of high consanguinity were rare and only observed among members of farming societies in our sample. Despite the lack of evidence for common consanguinity in our ancient sample, consanguineous traditions are today prevalent in various modern-day Eurasian societies,1,35-37 suggesting that such practices may have become widespread within the last few millennia.

Million-year-old DNA sheds light on the genomic history of mammoths.

Temporal genomic data hold great potential for studying evolutionary processes such as speciation. However, sampling across speciation events would, in many cases, require genomic time series that stretch well back into the Early Pleistocene subepoch. Although theoretical models suggest that DNA should survive on this timescale1, the oldest genomic data recovered so far are from a horse specimen dated to 780-560 thousand years ago2. Here we report the recovery of genome-wide data from three mammoth specimens dating to the Early and Middle Pleistocene subepochs, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these lineages gave rise to the woolly mammoth and the other represents a previously unrecognized lineage that was ancestral to the first mammoths to colonize North America. Our analyses reveal that the Columbian mammoth of North America traces its ancestry to a Middle Pleistocene hybridization between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were already present one million years ago. These findings highlight the potential of deep-time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Human population dynamics and Yersinia pestis in ancient northeast Asia.

We present genome-wide data from 40 individuals dating to c.16,900 to 550 years ago in northeast Asia. We describe hitherto unknown gene flow and admixture events in the region, revealing a complex population history. While populations east of Lake Baikal remained relatively stable from the Mesolithic to the Bronze Age, those from Yakutia and west of Lake Baikal witnessed major population transformations, from the Late Upper Paleolithic to the Neolithic, and during the Bronze Age, respectively. We further locate the Asian ancestors of Paleo-Inuits, using direct genetic evidence. Last, we report the most northeastern ancient occurrence of the plague-related bacterium, Yersinia pestis Our findings indicate the highly connected and dynamic nature of northeast Asia populations throughout the Holocene.

Temporal changes in the gene expression heterogeneity during brain development and aging.

Cells in largely non-mitotic tissues such as the brain are prone to stochastic (epi-)genetic alterations that may cause increased variability between cells and individuals over time. Although increased inter-individual heterogeneity in gene expression was previously reported, whether this process starts during development or if it is restricted to the aging period has not yet been studied. The regulatory dynamics and functional significance of putative aging-related heterogeneity are also unknown. Here we address these by a meta-analysis of 19 transcriptome datasets from three independent studies, covering diverse human brain regions. We observed a significant increase in inter-individual heterogeneity during aging (20 + years) compared to postnatal development (0 to 20 years). Increased heterogeneity during aging was consistent among different brain regions at the gene level and associated with lifespan regulation and neuronal functions. Overall, our results show that increased expression heterogeneity is a characteristic of aging human brain, and may influence aging-related changes in brain functions.

Molecular footprint of Medawar's mutation accumulation process in mammalian aging.

Medawar's mutation accumulation hypothesis explains aging by the declining force of natural selection with age: Slightly deleterious germline mutations expressed in old age can drift to fixation and thereby lead to aging-related phenotypes. Although widely cited, empirical evidence for this hypothesis has remained limited. Here, we test one of its predictions that genes relatively highly expressed in old adults should be under weaker purifying selection than genes relatively highly expressed in young adults. Combining 66 transcriptome datasets (including 16 tissues from five mammalian species) with sequence conservation estimates across mammals, here we report that the overall conservation level of expressed genes is lower at old age compared to young adulthood. This age-related decrease in transcriptome conservation (ADICT) is systematically observed in diverse mammalian tissues, including the brain, liver, lung, and artery, but not in others, most notably in the muscle and heart. Where observed, ADICT is driven partly by poorly conserved genes being up-regulated during aging. In general, the more often a gene is found up-regulated with age among tissues and species, the lower its evolutionary conservation. Poorly conserved and up-regulated genes have overlapping functional properties that include responses to age-associated tissue damage, such as apoptosis and inflammation. Meanwhile, these genes do not appear to be under positive selection. Hence, genes contributing to old age phenotypes are found to harbor an excess of slightly deleterious alleles, at least in certain tissues. This supports the notion that genetic drift shapes aging in multicellular organisms, consistent with Medawar's mutation accumulation hypothesis.

Ancient Mitochondrial Genomes Reveal the Absence of Maternal Kinship in the Burials of Çatalhöyük People and Their Genetic Affinities.

Çatalhöyük is one of the most widely recognized and extensively researched Neolithic settlements. The site has been used to discuss a wide range of aspects associated with the spread of the Neolithic lifestyle and the social organization of Neolithic societies. Here, we address both topics using newly generated mitochondrial genomes, obtained by direct sequencing and capture-based enrichment of genomic libraries, for a group of individuals buried under a cluster of neighboring houses from the classical layer of the site's occupation. Our data suggests a lack of maternal kinship between individuals interred under the floors of Çatalhöyük buildings. The findings could potentially be explained either by a high variability of maternal lineages within a larger kin group, or alternatively, an intentional selection of individuals for burial based on factors other than biological kinship. Our population analyses shows that Neolithic Central Anatolian groups, including Çatalhöyük, share the closest affinity with the population from the Marmara Region and are, in contrast, set further apart from the Levantine populations. Our findings support the hypothesis about the emergence and the direction of spread of the Neolithic within Anatolian Peninsula and beyond, emphasizing a significant role of Central Anatolia in this process.

Genomic and Strontium Isotope Variation Reveal Immigration Patterns in a Viking Age Town.

The impact of human mobility on the northern European urban populations during the Viking and Early Middle Ages and its repercussions in Scandinavia itself are still largely unexplored. Our study of the demographics in the final phase of the Viking era is the first comprehensive multidisciplinary investigation that includes genetics, isotopes, archaeology, and osteology on a larger scale. This early Christian dataset is particularly important as the earlier common pagan burial tradition during the Iron Age was cremation, hindering large-scale DNA analyses. We present genome-wide sequence data from 23 individuals from the 10th to 12th century Swedish town of Sigtuna. The data revealed high genetic diversity among the early urban residents. The observed variation exceeds the genetic diversity in distinct modern-day and Iron Age groups of central and northern Europe. Strontium isotope data suggest mixed local and non-local origin of the townspeople. Our results uncover the social system underlying the urbanization process of the Viking World of which mobility was an intricate part and was comparable between males and females. The inhabitants of Sigtuna were heterogeneous in their genetic affinities, probably reflecting both close and distant connections through an established network, confirming that early urbanization processes in northern Europe were driven by migration.